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ABSTRACT 

The accretion disk around a compact object is a nonlinear general rclativistic system involving mag- 
nctohydrodynamics. Naturally the question arises whether such a system is chaotic (deterministic) or 
stochastic (random) which might be related to the associated transport properties whose origin is still 
not confirmed. Earlier, the black hole system GRS 1915+105 was shown to be low dimensional chaos in 
certain temporal classes. However, so far such nonlinear phenomena have not been studied fairly well 
for neutron stars which are unique for their magnetosphere and kHz quasi-periodic oscillation (QPO). 
On the other hand, it was argued that the QPO is a result of nonlinear magnetohydrodynamic effects 
in accretion disks. If a neutron star exhibits chaotic signature, then what is the chaotic/correlation 
dimension? We analyze RXTE/PCA data of neutron stars Sco X-l and Cyg X-2, along with the black 
hole Cyg X-l and the unknown source Cyg X-3, and show that while Sco X-l and Cyg X-2 are low 
dimensional chaotic systems, Cyg X-l and Cyg X-3 are stochastic sources. Based on our analysis, we 
argue that Cyg X-3 may be a black hole. 

Subject headings: stars: neutron — X-rays: binaries — X-rays: individual (Sco X-l, Cyg X-l, 2, 3) — 
accretion, accretion disks 



1. INTRODUCTION 

X-ray binary systems vary on timescales rang- 
ing from mon t hs to milli-se c onds (see, e.g. , 
tiChen et al. 19971: iPaul et al. 19971 ; iNowak et al. 19991 : 
ICui 19991: IGleissner et al. 2004t lAxelsson 2"ool V. 
Detailed analysis of their temporal variability and 
fluctuation provides important insights into the geom- 
etry and physics of emitting regions and the accretion 
process. However, the origin of variability is still not 
clear. It could be due to varying external parameters, 
like the infalling mass accretion rate. It could also 
be due to possible instabilities in the inner regions of 
the accretion disk where the flow is expected to be 
nonlinear and turbulent. Uttley et al. (2005) (see also 
Timmer et al. 2000 and Thiel et al. 2001) argued that 
the non-linear behavior of a system can be understood 
from the log-normal distribution of the fluxes and 
the rms-flux relation. This implies that the temporal 
behavior of the system may be driven by underlying 
stochastic variations. By studying the underlying non- 
linear behavior, important constraints can be obtained 
on these various possibilities. 

An elegant way of obtaining the constraint is to per- 
form the nonlinear time series analysis of observed data 
and to compute the correlation dimension D2 in a non- 
subjective manner. Thi s technique has already been used 
to diverse situations (dGrassberger fc Procaccia 1 983a: 



Grassberger fe Procaccia 19 83b; 



Abcrbandcl 199 
Misra et al. 200 



Schrcibcr 199£ 



Serre et al. 19961 IMisra et al. 20041 ; 
Harikris hnan et al. 20061 ) . and refer- 
ences therein). By obtaining D2 as a function of the 
embedding dimension M, one can infer the origin of the 
variability. For example, D2 « M for all M corresponds 
to the system having stochastic fluctuation which favors 
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the idea that X-ray variations are driven by variations 
of some external parameters. On the other hand, a 
saturated D2 to a finite (low) value, beyond a certain 
M, implies a deterministic chaos which argues in favor 
of inner disk instability. However, to implement the 
algorithm successfully, the system in question should 
provide enough data. 

The technique was used earlier to understand the 
no nlinear nature o f a black hole system Cyg X- 
1 ju nno et al. 1990[) and an Active Galactic Nucleus 
(AGN) Ark 564 (jGliozzi et al. 200% . but due to in- 
sufficient data points the analyses were hampered 
and no concrete conclusions were made about £>2- 
Later on, anoth er black hole sys te m GRS 1915+105 
was analyzed (IMisra et al. 200l IMisra et al. 20061 ; 
lHarikrishnan et al. 20061 ) which was shown to display 
low dimensional chaos in certain temporal classes, while 
stochastic in other classes. 

However, so far none of the neutron star systems 
have been analyzed in detail in order to understand 
the origin of nonlinearity. Decades back, Voges et al. 
(1987) attempted to understand the chaotic nature of 
Her X-l, but th e analysis was hampered by low signal 
to noise ratio orris fc Matilskv 1989f ). Since then 
the investigation of chaotic signature in neutron stars 
remains unattended. Can a neutron star system not be 
deterministic? Indeed several features of X-ray binary 
systems consisting of a neutron star, such as their magne- 
tosphere and kHz Quasi-Periodic Oscillation (QPO) and 
its possible relation to the spin frequency of the neutron 
star, favor t he idea that they exhibit nonlinear res o- 
nance (e.g. (jBlaes et al. 20071: iMukhopadhvav 2009D 1. 
While the QPO itself is a mysterious feature whose 
origin is still unclear, its possible link to the spin 
frequency of the neutron star 4 indicates the origin 

4 However, some authors (Mendez & Belloni 2007) suggested 
that the kHz QPOs may not be related to the spin. 
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of QPO to be from nonlinear phenomena. Several 
LMXBs having a ne u tron star exhibit t win kHz QPOs 
(|Mendez et al. 1998t Ivan der Klis 20061). For som e 
of them, e.g. 4U 1636-53 dJonker et al. 2002Q . 
KS 17 31-260 (iSmith et al 19971) . 4U 1702- 
429 dMarkwardt et al. 19991) . 4U 1728-34 
(|van Straaten et al. 20021 ). the spin frequency of 
the neutron star has been predicted from observed data. 
However, for the source Sco X-l, which exhibits no - 
ticeable time variability (jMendez &: van der Klis 200(1 . 
while we observe twin kHz QP Os, we do not know t he 
spin frequency yet (but see (Mukhop adhvav 2 009)). 
For another neutron star Cyg X -2, we do observe 
kHz QPOs (jWiinands et al. 19981 ) as well. Several 
black holes also ex hibi t QPOs, e.g. GRS 1915+10 5 
dBelloni et al. 20011: iMcClintock fe Remillard 2006h . 
Cvg X-l dAngelini et al. 1994h . 

In the present paper, we first aim at analyzing the time 
series of two neutron star sources Sco X-l and Cyg X-2 to 
understand if a neutron star is a deterministic nonlinear 
(chaotic) system. Then we try to manifest the knowl- 
edge of nonlinear (chaotic/random) property of compact 
sources to distinguish a black hole from a neutron star. 
Subsequently, knowing their difference based on the said 
property, we try to identify the nature of a unknown 
source (whether it is a black hole or a neutron star). 
While the nature of some sources, as mentioned above, 
has already been predicted based on alternate method, 
for some others, e.g. Cyg X-3, SS433, it has not yet been 
confirmed. 

For the present purpose, we therefore concentrate on 
three additional sources Cyg X-l, Cyg X-2 and Cyg X-3. 
While Cyg X-l has been predicted to be a black hole and 
Cyg X-2 be a neutron star, t he nature of Cyg X-3 is not 
confirmed yet. Some authors (lErgma &: Yungelson 19981 : 
iSchmutz et al. 1996t iSzostek &; Zdziarski 20081 ) argued 
for a black hole nature of Cyg X-3, on the basis of its 
jet, the time variations in the infrared emission lines, the 
BeppoSAX X-ray spectra and so on. However, earlie r it 
was argued for a neutron star (jChadwick et al. 19851 ) by 
measuring its 1000 GeV 7-rays which suggests a pulsar 
period of 12.5908 ± 0.0003 ms. By analyzing the time 
series and computing the correlation dimension D 2 , here 
we aim at pinpointing the nature of Cyg X-3: whether a 
black hole or a neutron star. 

In the next section, we briefly outline the procedure 
to be followed in understanding the nonlinear nature of 
a compact object from observed data and to implement 
it to analyze the neutron star source Sco X-l. In §3, we 
then describe nonlinear behaviors of Cyg X-l, Cyg X-2 
and Cyg X-3. Subsequently, in §4, we compare all the 
results and argue for a black hole nature of Cyg X-3. 
Finally, we summarize in §5. 

2. PROCEDURE AND NONLINEAR NATURE OF SCO X-l 

Th e method to obtain D 2 is already established (see , 
e.g.. dGrassberger & Procaccia 1 983b: Mi sra et al. 20041 : 



lHarikrishnan et al. 2006h ). Therefore, here we discuss it 
briefly. We consider PCA data of the RXTE satellite 
(see Table 1 for the observations IDs) from the archive 
for our analysis. We process the data using the FTOOLS 
software. We extract a few continuous data streams of 
2500 — 3500 sec long. The time resolution used to gen- 
erate lightcurves is ~ 0.1 — 1 sec. This is the range of 



optimum resolution, at least for the sources we consider, 
to minimize noise without losing physical information of 
the sources. A finer time resolution would be Poisson 
noise dominated and a larger binning might give too few 
data points to derive physical parameters from it (see 
Misra et al. 2004, 2006, for details). 

Then we calculate the correlation dimension accord- 
ing to the Grassberger & Procaccia (1983a, b) algorithm. 
From the time series s(ti) (i = 1,2,.. .,N), we construct an 
M dimensional space (called embedding space ) , in which 
any vector has the following form: 

x(U) = [s(U), s{U + r), ,s{U + (M - l)r)], (1) 

where r is the time delay chosen in such a way that 
each component of the vector x(ti) is independent of each 
other. For a particular choice of embedding dimension 
M, we compute the correlation function: 



C M (r) 



N 



N c 



i=l j=1,j=£i 



\xi- Xj \), (2) 



which is basically the average number of points within a 
hypersphere of diameter r, where O is a Heaviside step 
function, N the total number of points and N c the num- 
ber of centers. If the system has a strange attractor, then 
one can show that for a small value of r 



D 2 (M) 



d log C M (r) 
d log r 



(3) 



In this numerical calculation, we divide the whole 
phase space into M cubes of length r around a point 
and we count the average number of data points in 
these cubes to calculate Cm(t). The edge effects, which 
come due to the finite number of data points, have been 
avoided by calculating C«(r) in the range r m i n < r < 
r-max, where r m i n is the value of r for Cm (r) just greater 
than one and r max can be found by restricting the M 
cubes to be within the embedding space. In Fig. [TJ we 
show the variation of log(Cjw(r)) with log(r) for different 
values of M for Sco X-l data. 

D2(M) can be calculated from the linear part of the 
log(Cftf(?*)) vs. log(r) curve and its value depends on the 
value of M. For a stochastic system, D 2 k, M for all M. 
On the other hand, for a chaotic or deterministic system, 
initially D 2 (M) increases linearly with the increase of M, 
then it reaches a certain value and saturates. This satu- 
rated value of D 2 is taken to be the correlation dimension 
of the system which is a non-integer. The standard de- 
viation gives the error in D 2 . 

We first concentrate upon the neutron star source 
Sco X-l. In Figs. [2K,b we show respectively the 
lightcurve and the variation of D 2 as a function of 
M and find that D 2 saturates to a value 2.6 (± 0.8). 
As this is a non-integer, the system might be chaotic. 
On the other hand, we know that the Lorenz attrac- 
tor is an example of an ideal chaos with D 2 — 2.05. 
Therefore, Sco X-l may be like a Lorenz system. But 
due to noise its D 2 seems appearin g slightly higher 
(|Misra et al. 20041 : iMisra et al. 2006D than the actual 
value. However, one should be cautious about the fact 
that Sco X-l is a bright source (much brighter than other 
sources considered later) . Hence, the dead time effect on 
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the detector might affect the actual value of saturated 
-D2 and the computed value might be slightly different 
than the actual one. However, this can not rule out the 
signature of chaos in Sco X-l, particularly because the 
corresponding count rates are confined in the same order 
of magnitude and hence the dead time effect, if any, is 
expected to affect all the count rates in a similar way. 

However, a saturated D2 is necessary but not a suf- 
ficient evidence for chaos. Existence of color noise 
(for which the power spectrum P{v) oc v~ a , where 
the power spectral indices a = 0, 1 and 2 correspond 
to "white" , "pink" and "red" noise respectively) into 
a stochastic system might lead to a saturated D2 of 
low value as well (e.g. Osborne & Provenzale 1989; 
Theiler et al 1992; Misra et al. 2006; Harikrish- 
nan et al. 2006). Therefore, it is customary to an- 
alyze data by alternate appro ach(s) to disting uish it 
from a pure noisy time series (jKugiumtzis 19 99). One 
of the techniques is the s urrogate data analysis (e.g. 
(jSchreiber fc Schmitz 19 96)). which has been described 
earlier in detail and implemented for a black hole 
(jMisra et al. 20061 ; IHarikrishnan et al. 200l . In brief, 
surrogate data is random data generated by taking the 
original signal and reprocessing it so that data has the 
same/similar Fourier power spectrum and autocorrela- 
tion along with the same distribution, mean and vari- 
ance as of the original data, but has lost all deterministic 
characters. Then the same analysis is carried out to the 
original data and the surrogate data to identify any dis- 
tinguishable feature(s) between them. The scheme pro- 
posed by Schreiber & Schmitz (1996), known as Iter- 
ative Amplitude- Adjusted Fourier Transform (IAAFT), 
is more consistent to generate surrogate data. 

Figures [2j:-f compare results for the original data with 
the surrogate data. It is clear that while distributions 
and power spectra are same/similar for both the data 
sets, Z?2 is much higher for the surrogate data which 
suggests existence of low dimensional chaos in Sco X-l 
with D2 ~ 2.6. This confirms, for the first time to best of 
our knowledge, a neutron star source to display chaotic 
behavior. As the existence of chaos is a plausible signa- 
ture of instability in the inner region of accretion flows 
which is nonlinear and turbulent, as mentioned in §1, the 
corresponding QPO, which is presumably an inner disk 
phenomenon as well, is expected to be governed by non- 
linear resonance mechanisms (e.g. Mukhopadhyay 2009). 

In Table 1 , we enlist the average counts < S >, its root 
mean square (rms) variation V < S 2 > — < S > 2 / < 
S >, the expected Poisson noise < PN > = V < S >, 
and the ratio of the expected Poisson noise to the rms 
value for all sources. It clearly shows a strong correla- 
tion between the inferred behavior of the systems and 
the ratio of the expected Poisson noise to the rms value. 

3. NONLINEARITY OF CYG X-l, 2,3 

We now look into three additional compact sources: 
Cyg X-l (black hole), Cyg X-2 (neutron star) and Cyg X- 
3 (nature is not confirmed yet), and apply the same 
analysis as in the case of Sco X-l. Figure [5b shows 
that Z?2 for Cyg X-l seems not to saturate and ap- 
pears very high. However, there is no surprise in it 
because its variability is similar to the temporal class 
X of the black hole GRS 1915+105 which was shown 
to be Poisson noise dominated and stochastic in nature 



(jMisra et al. 20 04). Indeed, earlier analysis of Cyg X-l 
data, while it could not conclusively quantify the under- 
lying chaotic behavior due to insufficient data, revealed 
very high dimensional chaos. Moreover a large < PN > 
(as well < PN > /rms) for Cyg X-l, compared to that 
for Sco X-l given in Table 1, reveals the system to be 
noise dominated. This ensures Cyg X-l to be different 
from Sco X-l. However, the variation of D2 as a func- 
tion of M for the original data does not deviate notice- 
ably from that of corresponding surrogate data, as shown 
in Fig. [3fc, which argues that Cyg X-l is not a chaotic 
system. 

Figures Hp, c show that D 2 for Cyg X-2 saturates to a 
low value ~ 4, which is significantly different than that 
of corresponding surrogate data. The power spectra and 
distributions, on the other hand, for original and surro- 
gate data are same/similar (as shown in Sco X-l, not 
repeated further). The saturated D2 for Cyg X-2 is al- 
most double than that of Lorenz system, possibly due to 
high Poisson noise to rms ratio (see Table-1). This sug- 
gests the corresponding system to be a low dimensional 
chaos. 

From Figs. 03, c we see that for Cyg X-3 the varia- 
tions of D2 as a function of M for original and surrogates 
data are similar to that of Cyg X-l. This confirms that 
the behavior of the unknown source Cyg X-3 is similar 
to that of the black h ole source Cyg X-l (see, however, 
(jAxelsson et al. 20081 )). Note from Table 1 that Cyg X-l, 
X-2, X-3 are significantly noise dominated compared to 
Sco X-l. Although noise could not suppress the chaotic 
signature in the neutron star Cyg X-2, its saturated D 2 
is higher than that of Sco X-l. On the other hand, even 
though the Poisson noise to rms ratio in Cyg X-l is lower 
than that in Cyg X-2 (but Poisson noise itself is higher 
in Cyg X-l), its D2 never saturates, which confirms the 
source to be non-chaotic; the apparent stochastic signa- 
ture is not due to Poisson noise present into the system. 

4. COMPARISON BETWEEN CYG X-l, CYG X-2 AND 
CYG X-3 

Finally, we compare the variations of D2 for all three 
cases of Cygnus in Fig. [6] Remarkably we find that 
D 2 values for Cyg X-l and Cyg X-3 practically overlap, 
appearing much larger compared to that for Cyg X-2 
which is shown to be a low dimensional chaotic source. 
On the other hand, Cyg X-2 is a confirmed neutron star 
and Cyg X-l a black hole. Therefore, Cyg X-3 may be a 
black hole. 

5. SUMMARY 

The source Cyg X-3, whose nature is not confirmed yet, 
seems to be a black hole based on the analysis of its non- 
linear behavior. On the other hand, we have shown, for 
the first time to best of our knowledge, that neutron star 
systems could be chaotic in nature. The signature of de- 
terministic chaos, which argues in favor of inner disk in- 
stability, into an accreting system has implications in un- 
derstanding its transp ort properties particu larly in Ke- 
plerian accretion disk (Win ters et al. 2003h . Note that 
in Keplerian accretion disks transport is necessarily due 
to turbulence in absence of significant molecular viscos- 
ity. The signature of chaos confirms instability and then 
plausible turbulence. On the other hand, for a rotating 
neutron star having a magnetosphere, signature of chaos 
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sugg ests their QPOs to be nonlinear resonance phenom- 
ena (|Mukhopadhvav 2009f) . The absence of chaos and 
related/plausible signature of instability in Cyg X-l and 
Cyg X-3 sugge sts the underlying a c cretion disk to be 
sub-Keplerian (jNaravan & Yi 19951 : iChakrabarti 1996f ) 
in nature which is dominated significantly by gravita- 
tional force. 
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Fig. 1. — Variation of log (Cm) as a function of log( r) for different embedding dimensions. The linear scaling range is used to calculate 
the correlation dimension. 

TABLE 1 

Observed data 



Source 


Obs. I. D. 


< S > 


rms 


< PN > 


< PN > /rms 


Behavior 


Sco X-l 


91012-01-02-00 


58226 


0.074 


0.004 


0.054 


C 


Cyg X-l 


10512-01-09-01 


10176 


0.261 


0.031 


0.119 


NC/S 


Cyg X-2 


10063-10-01-00 


4779 


0.075 


0.014 


0.191 


C 


Cyg X-3 


40061-01-07-00 


3075 


0.125 


0.057 


0.455 


NC/S 



Columns:- 1: Name of the source, 2: RXTE Observational identification number from which the data has been extracted. 3: 
The average count in the lightcurve < S > 4: The root mean square variation in the lightcurve, rms. 5: The expected Poisson noise 
variation, < PN >= V< S >. 6: The ratio of the expected Poisson noise to the actual root mean square variation 7: The behavior of the 
system (C: chaotic behavior; S: stochastic behavior; NC: nonchaotic behavior) 
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Fig. 2. — Sco X-l: (a) Variation of count rate as a function of time in units of 10 7 sec (lightcurve) , without subtracting the initial 
observation time, (b) Variation of correlation dimension, along with error bars, as a function of embedding dimension for original data. 
The solid line along the diagonal of the figure indicates an ideal stochastic curve, (c) Variation of correlation dimension as a function of 
embedding dimension for original (points) and corresponding surrogate (dashed lines) data, (d) Variation of number of count rate as a 
function of count rate itself in units of 10 4 sec -1 (Distribution) for original (solid line) and surrogate (points) data. Power-spectra for (e) 
original and (f) surrogate data. 
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Fig. 4. — Cyg X-2: (a) Variation of count rate as a function of time in units of 10 7 sec (lightcurve). (b) Variation of correlation dimension, 
along with error bars, as a function of embedding dimension for original data. The solid line along the diagonal of the figure indicates an 
ideal stochastic curve, (c) Variation of correlation dimension as a function of embedding dimension for original (points) and corresponding 
surrogate (dashed lines) data. 




Fig. 5.— 



Cyg X-3: Same as Fig. \4\ 



